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Abstract 

We study the transition form factor of 7r° — > 7*7 as a function of the momentum transfer Q 2 
within the light-front quark model (LFQM). We compare our result with the experimental data 
by BaBar as well as other calculations based on the LFQM in the literature. We show that our 
predicted form factor fits well with the experimental data, particularly those at the large Q 2 region. 
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I. INTRODUCTION 



The BaBar collaboration [l| has reported a new data of the 7r° — > 7*7 transition form 
factor F nj (Q 2 ) for the high momentum transfer Q 2 up to 40 GeV 2 . To describe the data 
with the Q 2 dependence, the form factor is fitted to satisfy the formula 



Q 2 \F^{Q 2 )\=A 



Q 2 



10 GeV^ 



(1) 



with A = 0.182 ± 0.002 GeV and /3 = 0.25 ± 0.02. Before the new data, most theoretical 
models predicted that the form factor approaches the QCD asymptotic limit 2J, depending 
on the pion distribution amplitude (DA) with the Q 2 dependence under 10 GeV 2 3j-|5|. 
Obviously, the experimental values for Q 2 > 10 GeV 2 by BaBar are surprisingly much 
higher than the QCD asymptotic expectations and thus, cannot be explained by the lowest 
perturbative results Even the high order corrections are considered BE], the large Q 2 
behavior is still hard to be understood. Recently, many proposals js 31j have been given 
in the literature to understand the transition form factor, particularly the BaBar data for 
Q 2 > 10 GeV 2 . 

In this note, we will use the phenomenological light front (LF) p ion wave function to 
evaluate Q 2 \F ni (Q 2 )\ in the light front quark model (LFQM) 32j-|36[]. We will concentrate 
on the space-like region for the transition form factor. The LF wave function is manifestly 
boost invariant as it is expressed in terms of the longitudinal momentum fraction and relative 
transverse momentum variables. The parameter in the hadronic wave function is determined 
from other information and the meson state of the definite spins can be constructed by the 
Melosh transformation. We emphasize that our derivation of the form factor can be applied 
to all allowed kinematic region. In Ref. 371] , the study on the transition pion form factor 
based on the LFQM has been done but the calculation for Q 2 is only up to 8 GeV 2 . With 
the same set of parameters in Ref. [3?]], the high Q 2 BaBar data cannot be fitted. The use 
of the LFQM to understand the BaBar data has been explored in Ref. {3^ . However, the 



conclusion in Ref. [38( has failed to explain the data. In this work, we would like to revisit 
the LFQM to see if it is indeed the case. 

This paper is organized as follows. In Sec. II, we present the relevant formulas for the 
matrix element and form factor for the 7r° — > 7*7 transition. In Sec. Ill, we show our 
numerical analysis. We give our conclusions in Sec. IV. 



II. THE FORM FACTOR 



The transition form factor of F 7r o_ >7 * 7 «(g , f , which describes the vertex of 7r°7*7*, is 
defined by: 

A(tt°(P) -+7*(<Zi,ei) 7*fe,e 2 )) = ie 2 F^ rr (qlq 2 2 ) W ^ <tf g 2 CT , (2) 

where F 7r o_ 5 . 7 * 7 *(g^, gf) is a symmetric function under the interchange of q\ and q\. From the 
quark- meson diagram depicted in Fig. 1, the amplitude in Eq. <^ is found to be 




0(P) P3,u 7 *(g 2 ) ^(p) Pi," 7 * fe) 

(a) (b) 



FIG. 1. Loop diagrams that contribute to n c 
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(3) 



where N c is the number of colors, cq is the quark electric charge and Ap is the vertex function 
related to the 7r° meson bound state. To calculate the 7r° — > 7*7* transition from factor 
within the LFQM, we have to decompose the 7r° meson into QQ Fock states, described as 
(wh — dd) I \/2. In the LF approach, the LF meson wave function can be expressed by an 
anti-quark Q and a quark Q with the total momentum P as: 



\M(P, S,S Z )) = Y.J [<fod [dp2}2(27T) 3 6 3 (P - Pl - pa) 
A1A2 



x$^(z,fc ± )6 ( t( Pl ,A 1 )4(p 2 ,A 2 )|0) 



where 



[dp] 



dp + d 2 p± 
2(2vr) 3 



(4) 



(5) 



$^ A2 is the amplitude of the corresponding q(q) and p\M) is the on-mass shell LF momentum 



of the internal quark. In the momentum space, the wave function $^f z is given by 



$ s J'(h, k 2 , Ai, A 2 ) = RfX(. z > fc -0 <f>& M 



ss z 



(6) 



where <f)(z, k±) represents the momentum distribution amplitude of the constituents in the 
bound state and R\ \ 2 constructs a spin state (S, S z ) out of light front helicity eigenstates 



(AiA 2 ) 



The LF relative momentum variables (z, k±) are defined by 
p+ = Z P + , 



p+ = (1 - z)P^ 



Pi± = zP± - k±, p 2 ± = (1 - z)P± + k ± . 

The normalization condition of the meson state is given by 

(M(P', S', S' Z )\M(P, S, S z )) = 2(2vr) 3 P+5 3 (P f - P)6 s , s 6 s , Sm 

which leads the momentum distribution amplitude <j>(z, k±) to 

dz d 2 k 



i . 



2(2vr) 3 

We note that Eq. (JSJ) can, in fact, be expressed as a covariant form 



(7) 



(8) 



(9) 



32 



33 



40] 



qSS z 



M 



vivi 



Ml 



2[M*-(m Q -m Q y 
m Q + k l , rn 2 Q + k 2 



u (Pi> ^i) 7 5 ^ (P2, A 2 ) 4>(z, k±) 



(10) 



z 1 — z 

In principle, the momentum distribution amplitude <j>(z, k±) can be obtained by solving 



the LF QCD bound state equation 



331 ] . However, before such first-principle solutions are 



available, we would have to be contented with phenomeno 



that has been used is the Gaussian type wave function 

Mz,k\ 



ogical amplitudes. One example 



3J-[36|: 



1 dk z 
A/ \l — — exp 



N r dz 



2ul 



(11) 



where N = 4(ir / cu^) ± , k = (k±, k z ), and k z is defined through 

En + k z 



En + Er 



1-z 



Eq — k z 

En + En 



E; 



mf + k 2 



(12) 



by 



and dk z /dz = EqEq/z(1 — z)M . After integrating over p 3 in Eq. (J3]), we obtain 



A(QQ -> 7*(?0 7*fe)) = e Q e Q N c j q J dpj J 



d 2 p 3 ± 



A, 



12 -P2c 



+ (ei e 2 , gx g 2 ) 



+ (Pl(3) ^ P3(l)) , (14) 



and 



/ = Tr[7 5 (- fa + mg) $(#2 + mg) + toq)] , 



Pi 



■ + Pa 
Pi 



(15) 



where the subscript {on} stands for the on-shell particles. One can extract the vertex 



function Ap from Eqs. (J3J), f fTOj) and (Jl4j) . given by 

Ap 



32 



Pi Pa 



(2;, fc ± ) 



(16) 



To calculate the trace /, we use the definitions of the LF momentum variables (z(x), kj_(k' ± )) 
and take the frame with the transverse momentum (P — q 2 )± = for the QQ state (P) and 
photon (q 2 ) in Fig. la. Hence, the relevant quark variables are: 

P t = zP + , P + = (1-Z)P + , Pl± = zP ± -k ± , P 3± = (l-z)P ± + k ± . 

Pt = xq£, pt = (1 - x)q£, p 2 ± = xq 2± - k' ± , p 3± = (1 - x)q 2± _ + k' ± . (17) 



At the quark loop, it requires that 



[z - x)q 2± + k\ 



(18) 



Take the trace I into Eq. ( jl~4l) and consider the 7r° meson Fock states, the form factor 
F w o_>. 7 . 7 * (qf, q 2 ) in Eq. (J2J) can be found to be: 



F 0^ * 

1 7T U — >7 7' 



/ 2 2^ 
(9l>92 



4 /A- f dx d 2 k 




6 



2(2 



7T 



, , 2 \ m Q + (1 - z)ro Q fcje 
l ' lj 4l-zk 2 - (m 2 + k 2 



+ (g 2 ^ gi) ^ + (Q 0) , 



(19) 



with 



' z{ 1 — z) dk 



dz 



2M 2 



exp 



2u M 



k — (k±,k z ), x = zr , 6 



A; 2 



$(z,fc 2 ; 



dki 



^ + 9i 2 ) + V( m l + Q2- 9i) 2 - 4( ?fraf 
2wip 



(20) 
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III. NUMERICAL RESULT 



To numerically evaluate the transition form factor of ir° — > 7*7, we need to specify the 
parameters in Eq. f[T9~j) . To constrain the quark masses of m u ^ s and the pion scale parameter 
of u n , we use the meson decay constant f w o and the decay branching ratio of 7r° — > 27, given 
by the PDG [42 1 



fa = 130 MeV, B(n° ->• 2 7 ) = (98.832 ± 0.034) 



where the explicit expressions of f w o 



fa 



43| and B(n° -» 2 7 ) are 
dx d 2 k 



V2 J 2(2tt) 



(x, k 1 



to 



^Jm 2 + k\ 



and 



-> 2 7 ) = ^?^|F(0,0)^ 27 | 2 



(21) 



(22) 



(23) 



64^ 

respectively. As an illustration, we extracte \F(0,0) n o^ 2 "/\ = 0.274 in GeV -1 , to = m u = 
ma = 0.24 and u n = 0.31 in GeV, which will be used in our following numerical calculations. 

We now consider the case with one of the photons on the mass shell. From Eq. f|T9|) . the 
transition pion form factor becomes 

4\/2 [N~ c [ fdxd 2 k_ L . l9N m + (1 - z)mk 2 6 



F-n-~f{Q ) — F w o^y*^(Q , 0) 



3 V 3 



2(2vr) 



TO + (1 



z{l - z)Q 2 - (to 2 + k\ 
z)mk 2 l G 



(24) 



2(2tt) j v ' x ' K + A; 2 ) 
In Fig. |2l we show the form factor in Eq. fl24|) . We note that the first term in Eq. (J2! 
dominates for the lower re gion of Q 2 and thus, it can be use to describe the experimental 
data of BaBar [1], CLEO 44] and CELLO 45^ with Q 2 < 10 GeV 2 . The second one in 
Eq. ( 1241) . related to the non- valence quark contributions, is small for a small Q 2 , but it may 
enhance the form factor with a high value of Q 2 . As a result, we will include this term in 
our our numerical calculations. To easily examine the Q 2 dependence of the form factor, we 
have fitted our result in terms of the double-pole form: 



F n ry (Q 



7T— »77 



(0,0) 



M+(/3Q) 2 -(aQ) 4 "" (25) 
Explicitly, we find that the dimension parameters of a = 0.325, = 1.15 and ^0^.^(0, 0) = 
0.274 in GeV -1 , and the dimensionless parameter of M = 3.6. In Fig. [3j we concentrate on 




FIG. 2. (Color online) Q 2 dependence of F 7r7 (Q 2 ) in the LFQM. 



the behavior of the form factor in the region with Q 2 < 10 GeV 2 . It is easy to see that our 
results fit the data well in this region similar to other theoretical calculations as expected. 
In Fig. HI we show the DA, <f)(z), as the function of the momentum fraction of the internal 
quark and meson longitudinal momenta, z, obtained by the integration of k± in Eq. (II ip . 

As we mentioned in the introduction, the BaBar result cannot be fitted by extending the 



study in Ref. 



371 ] to a high Q 2 . The main reason is due to the choices of the free parameters, 



such as the quark masses and u n , leading to a different sharp of the pion wave function. 



comes from DAs. In 



Similarly, the main difference between our results and those in Ref. 
particular, our DA shown in Fig. H] appears to be much broader. We note that in Refs. |8Nl2|. 
a more broader DA of the pion is utilized to fit the BaBar data, particularly in the high Q 2 
region. The results in these models differ slightly from ours only at large values of Q 2 . It is 



interesting to point out that our result is almost identical with that in the Regge model 21 1 



28 



Finally, we remark 



and the double logarithmic behavior from the chiral anomaly effects 
that the single data point at Q 2 = 27.31 GeV 2 by BaBar, which is apparently consistent 
with the QCD asymptotic limit, cannot be explained by this work within the framework of 
the LFQM. 




FIG. 3. (Color online) F 1T1 (Q 2 ) for Q 2 < 10 GeV 2 in the LFQM. 




FIG. 4. 4>(z) as a function of z in the LFQM. 

IV. CONCLUSIONS 

We have studied the form factors of tt° — > 7*7 within the LFQM. In our calculation, 
we have adopted the Gaussian-type wave function and evaluated the form factors for the 
momentum dependences in the all allowed Q 2 region. We have also parametrized the form 

8 



factor in terms of the double-pole form. Our numerical values are close to the experimental 
results by BaBar. In particular, our results of the transition form factor fit well with the 
experimental data in the high Q 2 region, which cannot be explained in the previous attempts 
based on the framework of the LFQM. Finally, we remark that due to the large uncertainty 
in the high Q 2 region for the BaBar data, further theoretical studies as well as more precise 
experimental data are clearly needed. If some future experiment could not confirm the 
BaBar data but be rather in agreement with the QCD asymptotic limit, the parameters of 
the LFQM in this study should be either modified or ruled out. 
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